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t-h ! Abstract 

> ; 

^ ■ We study a dissipative gauge theory of nonrelativistic fermions in 2+1 dimensions at 

ps) ■ zero temperature by the Wilsonian renormalization-group method. In this theory, we 

\Q ' incorporate in the fermion propagator a new term of the form in ■ sgn(u;) (where k is a 

0\ : 

parameter and u is the fermion frequency), which is usually induced by impurity effects. 

a ; 

In the previous papers, we studied this system for k = 0, and showed that there exists a 
non- Fermi-liquid infrared fixed point. In this paper, we address the question whether this 
non- Fermi-liquid behavior remains stable or not in the presence of impurity effects, i.e., 
the k term. Our results show that the non-Fermi-liquid fixed point is unstable for k^O 
^ ■ and an effective gauge coupling constant tends to vanish at low energies. However, in 

h : 

intermediate energy scales, the behavior of correlation functions for k ^ is controlled by 
the non- Fermi-liquid fixed point at k — 0, i.e., a crossover phenomenon appears. Physical 
implications of this phenomenon are discussed. 
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1 Introduction 



In condensed matter physics, low-dimensional electron systems are of current interest, 
which include the high-T c superconductivity, the fractional quantum Hall effect, and the 
recently discovered new spin ladder systems which exhibit superconductivity. Theory 
of nonrelativistic fermions interacting with a gauge field plays an important role as an 
effective reference theory for these interesting systems. A typical example of such gauge 
theory is found in the work of Halperin, Lee, and ReadQ, in which they studied a system 
of electrons in the half-filled Landau level. This gauge theory is also regarded as an 
important constituent to compose a low-energy effective field theory for the t-J model@. 

Recently, it has been recognized that the renormalization-group (RG) method can be 
very efficiently applicable to a system of nonrelativistic fermions for investigating its low- 
energy behavior 0. Nayak and WilczekQ] made a RG analysis of a kind of e-expansion 
for the above gauge theory in Ref. at zero temperature, and claimed existence of a 
nontrivial infrared (IR) fixed point in some parameter region. Following their works there 
appeared some related works |J. This fixed point is very interesting because it implies 
a departure from the conventional Fermi liquid theory; At that fixed point, the theory 
exhibits non-Fermi-liquid-like behavior similar to that of a marginal Fermi liquid for the 
high-T c superconductivity]?! or of a Luttinger liquid in one dimension ||. 

In this theory, a dissipative term is generated in the gauge-field propagator via quantum 
fluctuations of fermions, and this term becomes important at low energies. Thus its 
effect to the above nontrivial fixed point should be clarified. In the Nayak and Wilczek 
approach, the relevance of this term to the above fixed point is not reflected well, since 
the dissipative term has no divergence in their e-expansion, hence not being treated as 
a relevant coupling. In the previous paper 0, we considered a two-dimensional model of 
nonrelativistic fermions interacting with a gauge field having a dissipative term from the 
beginning, and carried out the Wilsonian RG (WRG) analysis, treating this dissipative 
term as a genuine running couping constant. We found that there still exists a nontrivial 
fixed point. In this sense, the fixed point found by Nayak and Wilczek is stable against 
the dissipative term. 
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Since this fixed point is so important, it is worth studying its stability from a more 
general point of view. In this paper, we study its stability against "impurities" in detail. 
Explicitly, we incorporate an additional term in the fermion propagator of the previous 
dissipative gauge theory of Ref. 0, and study its low-energy effects by the WRG method. 
The term added is 

—in ■ sgn(o;), 

where k is a parameter, u is the fermion frequency, and 

f 1 > 0) 

sgn(x) = 1 (1.1) 
[-1 (x < 0) 

Motivation for adding the above term is twofold. First, this term is just generated via 
scattering of fermions by impurities. Real physical systems necessarily contain impurities, 
and they often play an important role; for example, in localization phenomenon in the 
quantum Hall effect. It is natural to extend the previous analysis in Ref.JJ to a system 
containing impurities. Usually, one treats effects of random impurities by letting them to 
interact with fermions through a potential v (x), so the interaction Hamiltonian is given 
by 



H int = J dx^(x)Y^v(x - Xi)?p(x), (1.2) 



where ip(x) is the fermion field operator and X{ are coordinates of impurities. Practical 
calculations at one-loop level show that the self energy Y>k(uj) of fermions has the following 
signature term 

Efc(ic<j) = —in ■ sgn(cj), (1.3) 

where k(> 0) is related with the Fourier modes of v (x). In the classical equation of motion 
of a fermion, the resistive (friction) force is given by oc — kv, where v is the velocity of 
fermion. In the random phase approximation (RPA) with respect to impurities, the 
fermion propagator is simply given as 

G (k, uj) = — — ^— -, (1.4) 

too + tKsgn{uj) + fi — hk 

where is the energy of fermion with momentum k, and /i is the chemical potential. 
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The constant k(> 0) has another physical significance. This term modifies the mo- 
mentum distribution of the fermion density which exhibits a step-function behavior at 
zero temperature T = and k = 0. Actually, from the propagator ( |1.4j ), the density of 
fermions N(Ek) at E^ is obtained as 

r°° du 1 1 

k J-oo 2ir iu) + z/«sgn(a;) + // — ^ 2 

= sgn(ii-E fc )(i - -arctanr ^— — )+-. (1.5) 



The curve N(Ek) is plotted in Fig.l. From (|1.5| ) and Fig.l, we see that the k term removes 
the singularity at E^ — 0. 

As is well known, this smooth momentum distribution ( p..5|) looks similar to that of a 
system of fermions at finite temperature (T). The second motivation is related with this 
point. When we consider thermal effects at finite T in this model, quantum corrections 
induce this signature term in the fermion self energy, even if we start with k = 0|T0| . 
Thus we must consider this signature term from the beginning when we investigate the 
system at T > by using the WRG analysis, treating k as a running coupling constant. 
These motivations let us start with the fermion action with the signature term as given 
by Eq-O in Sec.2. 

As we explained above, in the case of k = there is the nontrivial IR fixed point 
at which the system exhibits non- Fermi-liquid-like behavior || |J. We want to study 
the stability of this non-Fermi-liquid fixed point and the RG flow in the presence of the 
signature term. Throughout the paper we consider the case of zero temperature (T = 0). 



The other important case of finite T will be discussed separately [10|. 

This paper is organized as follows. In Sec.2, we shall introduce the model. In Sec. 3, the 
RG prescription is explained. We shall employ some specific momentum regularization, 
which is suitable for the present gauge-fermion system. In Sec. 4, we shall calculate one- 
loop corrections of fermions and gauge bosons by the momentum-shell integration. In 
Sec. 5, renormalization constants and RG equations are explicitly obtained. Qualitative 
and numerical solutions of the RG equations are given. It is shown that the non-Fermi- 
liquid fixed point at k = is unstable for nonzero k. However, at certain intermediate 
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energy scales, the behavior of the system is controlled by that fixed point. This means 
that the system exhibits a cross-over phenomenon. In Sec. 6, by solving the RG equation 
for the fermion propagator, we obtain its explicit form at intermediate energy scales and 
also in the low-energy limit, confirming the above conclusion of the cross-over. Section 7 
is devoted for concluding remarks. 



2 The Model 

In this section we explain our model in detail, which describes a 2-dimensional system of 
nonrelativistic spinless fermions ip(x, r) interacting with a dissipative gauge field Ai(x, r) 
(i = 1,2), where < r < oo is the imaginary time. Relationships of this model with 
electrons in the half-filled Landau level and with strongly correlated electron systems for 
high-T c superconductivity are discussed in Refs.[p], Qj. 

To study the low-energy excitations around the Fermi surface and perform the WRG 
analysis for the gauge-fermion system, it is useful to introduce a momentum cutoff as 
shown in Figs. 2 and 3@. The region around the Fermi surface is divided into N segments, 
and the fermion field in each divided part is labeled by a, a = 1, N, as ip(a; p, to). The 
fermion momentum p is measured from a point on the Fermi surface k Fta , the center of 
the a-th part. The momentum cutoff A is introduced as \p\ = \k — k F>a \ < A, where k is 
the original fermion momentum. 

In the standard path integral formalism, the partition function of this model is given 

by 

Z = J[d$][diP][dA]exp(-S^-S A -S int3 ). (2.1) 

In ( |2.1| ), is the free part of fermions given in the Fourier (momentum and energy) 
representation by 

Sjf, = dp "^ijj(a;p,u)(—iLU — iKsgn(u) + VFe a ■ p)vjj(a-,p,uj), (2.2) 

J J a 

where e a is the unit vector in the direction of k F ^ a and 

d 2 p du 
"P=77T~^' du = —, v F = k F m. (2.3) 



The free part of gauge field Sa is given by 

Sa = [de[dqA{-q,-e)A b - 1 v B (q 2 - b + A 2 - b X-^-)A{q,e). (2.4) 

The second term, with the coefficient A > 0, is the so-called dissipative term (Landau 
dumping factor) that is induced by the vacuum polarizations of fermions at one-loop level. 
It is shown in the following discussion that the parameter A behaves nontrivially under 
RG transformation. We have taken the Coulomb gauge d{Ai = 0, then the gauge field Ai 
has only one dynamical degrees of freedom, A(q, e) = ie^qiA^q, e)/q. We have introduced 
a parameter 6 which controls fluctuations of the gauge field. This generalization was first 
introduced in Ref . JTJ . The case 6=1 corresponds to the Coulomb potential for electrons 
in the half-filled Landau level, and the case 6 = corresponds to the case of the t-J model. 
The three-point interaction term of the action Si n t$ is given by 

Sint,3 = igvF / dude l dp dq V — A(q, e)ip(a; p + q, u + e)ip(a; p, u), (2.5) 

J J a Q 

where p x q = e^piq^. In a nonrelativistic fermion-gauge field system, the action also 
contains the four-point interaction term like AAtpip. We shall neglect this term, because 
this term contributes only to the renormalization of the chemical potential of the fermion 
and the effective mass of the gauge field. Actually, in the field theoretical approach |5|], 
it is shown that the mass of gauge field vanishes because the contribution from the four- 
point interaction cancels that from the three-point interaction, as guaranteed by the gauge 
invariance. In the following discussion, we shall not discuss the contribution from these 
terms to the renormalization of the chemical potential and the gauge field mass. 

3 Wilsonian renormalization group program 

In this section, we shall explain the WRG program for our gauge-fermion system intro- 
duced in Sec. 2. In the previous paper 0, we studied the case k = by the WRG and 
showed that there is a nontrivial IR fixed point for 6 < 1. Let us explain the procedures 
of WRG program briefly once gain. It consists of the following three steps (i)-(iii). 

Step (i) Integrate over the high-momentum modes (modes with the momentum A/e* < 
p < A; £(> 0) is a scaling parameter) of fermion and gauge fileds. Then the action of the 
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low-momentum modes is obtained as follows; 

f M/e* 

= / du / dp 



x ^2ip(a;p,uj)R(t)(-iuj - iR K (t)/tsgn(^>) + R VF (t)v F e a ■ p)ip(a;p,oj) 

a 

S A = JdeJ^ dqA(-q,-e)A b - 1 v B (q 2 - b + A 2 - b R x (t)\^-)A(q,e), 

f fA/e* 
Sint,3 = / / 



A/e* 

e a x q 



x H iR g (t)gv F — A(q, e)ip(a; p + q,uj + e)ip(a; p, u), (3.1) 

a Q 

where R(t), etc. are t-dependent renormalization constants. We have neglected terms 
that are higher -orders in momentum or fields as irrelevant terms. This point is discussed 
somewhat in detail in Ref . fl 1[| . 



Step (ii) The original action S and the new action S = + Sa + S intt 3 are defined 
in the two different kinematical regions, and then we rescale momenta in order to return 
the new momentum cutoff A/e* to the original value A. Similarly, we require so that 
certain quadratic terms in the effective action have the same coefficients with the original 
ones by rescaling fields. At the same time we scale energy, momentum and the fields. In 
nonrelativistic theories like the present model, there is an ambiguity in assigning scaling 
dimensions of energy and fields, because there is no principle for that. This is in sharp 
contrast with relativistic systems where it is natural to assign same scalings for energy and 
momentum. Therefore we introduce a general scaling law which contains a free parameter 
£ > as follows, 

uj — >• uj = e^uj, 

p _ p = e * P , ( |p| < A/e*, \p\ < A) 
V>(a; p, uj) -» 4>{a\ p, uj) = e~^ +1)t ^R(t)*jj(a; p, uj), 
A{p,uj) -> A{p,uj) = e- {2+ ^ )t A{p,uj). (3.2) 

Step (iii) The effective action S is expressed in terms of the above variables; 

S^, = JduJ dp^^(a\p,u)(-iuj - i«(t)sgn(o)) +v F (t)e a ■ p)-ijj(a;p,uj), 



Sa = Jde J A dqA(-q, -e)A b - 1 v B (q 2 - b + A^ b \(t)-^)A(q, e), 
S in t,3 = J dude J dpdq^2ig(t)v F (t)-^-: — A(q, e)ip(a; p + q, Co + e)ip(a;p, u)).(3.3) 



rip d<? y^jgi^vpit) — 

Q 

The effective parameters (coupling constants) which flow as the scaling parameter t varies 
have been defined as follows, 



a(t) = e^)* Ra{t) a 

v F (t) = eWRvrftvp, 
A(t) = e^R x (t)\, 

K {t) = e^R K (t)K. (3.4) 

The simplest assignment of scaling dimensions is £ = I. One- loop calculation shows 
that RG recursion equations of the effective parameters depend on the values of £. But 
this ^-dependence must be superficial. Actually, as we show in later discussion, physical 
quantities are independent of £. 

4 One-loop quantum corrections and renormalization constants 

In this section, we shall perform the momentum-shell integration for an infinitesimal RG 
transformation with exp(rft) ~ 1 + dt at one- loop level and calculate the renormalization 
constants R, R Vp , R g , R\ and R K . With this result of infinitesimal RG transformation, we 
derive the RG equations of our running coupling constants, g(t),VF(t), X(t) and n(t). We 
consider the fermion part in the action in Sec. 4.1, the gauge-boson part in Sec. 4. 2, and 
the vertex part in Sec. 4. 3. 

4.1 Renormalization of the fermion part 

The bare action is renormalized by the one-loop quantum correction which is obtained 
by integrating over the high-momentum modes of fermion field and gauge field in the 
momentum region A(l — dt) < p < A. The effective action is obtained as follows 

— / cko / dpip(a; p, uo) —ito — msgn(u;) + VFe a ■ p + E(a; p, to) ^p(a; p, a;), (4.1) 
„ J jo 



where E(a;p, 10) is the fermion self-energy coming from the high-momentum modes. It is 
given by 

de / dg — £>(,(<?, e)G (a;p + g,u; + e), (4.2) 

JA(i-dt) — \ q J 



where D (q, e) is propagator of the gauge field obtained from ( [O] 

1 



D (q,e 



i(q)(\e\ + F( q )y 

( \ \ VbA 
7(9) = A , 

= 7^(g)A 6 -W- fe = ^(|) 2 - fe . (4-3) 

Similarly, G (a;p, is the fermion propagator, obtained from of fl2T2| ) as 

G (a;p,^) : 



icu + msgn(u;) — F p ' 

E p = v F e a ■ p, (4.4) 

with k > 0. 

The momentum-shell integration over the radial direction in Eq.( |4.2j) can be easily 
performed to obtain the self-energy as 

S(a;p,c<j) = —iaF^dt J dy\ sin (f)(y)\I (uJ, E,k), (4.5) 

where 



00 



I(uJ, E,1Z)= de ._ — . — • . (4.6) 

J -00 e + io + Ksgn(e + 10) + lE |e| + 1 



In the above we have introduced 



a = 



2n 2 v B ' 



— w F e a • g w F A cos ( 

E = y + — - — , y = 



e = 



F A ' " F A ' 

— , uJ= — , k = — . (4.7) 

-TA -TA -TA 



The (p is the angle between the two vectors p and e a . In the above momentum integration, 
we changed the integration variable from <f> to y. 
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We divide the integrand / ( |4.6| ) into six terms as there are two signature terms sgn(e+u;) 
and lei. 



3 3 

I(U,E,K) = 0(5?) I%(P,E,k) — &(—uj) Ii(—u3, -E,K) 



i=i 



i=i 



where 



I 2 (uj,E,k) 
I 3 (uJ,E,k) 
and 0(x) is the step function, 



de- 



' e + lu — n + iE —e+1' 







de- 



e + uJ + K + iE -e + 1 

1 1 
de- — — — — = • = , 

e + u + k + iE e+1 



Q(x) 



r 1 (x > 0) 

I (x < 0) 
The integrations in (|4.9j) can be evaluated as 



(4.9) 



(4.10) 



h(uj,E, k) 
I 2 (uj,E,k) 



1 



1 ln (ZJ+1] 



-7T + isgn(E)( arctan -= 

(ZJ+k) 2 + e" 



z 1 + 1 



-(In 

2 l k 2 + E 



— 2 + ln |u7 + l\ 2 ) 



— ijJ + K K 

-isgn{E) (arctan _ arctan -=- 



h(u,E, k) 



zi-1 



2 _ 2n , / _ N ,7T ; + K x 



-(ln((tJ + 7t) 2 + £ ) + isga{E)(- - arctan — = 
2 2 \hi 



where 



Z\ = u -\- iE + k 1 z 2 = uj + iE — k. 



(4.11) 



(4.12) 



We are interested in the renormalization constants, so we expand the self-energy 
E(a; p, lo) in powers of ioj and E p = VFt a ■ p and keep only the linear terms; 



S(a; p, = S + S w • iu + S p • E p . 



(4.13) 



Detailed calculations are given in Appendix A, and the expansion coefficients are obtained 
as follows; 



0. 
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S w = —adtKi(K v ,X), 

S p = -adt(K 2 {K v ,\) -K^X)), (4.14) 
where we introduced the parameter k v , a rescaled k, and two functions Ki t2 (n v , A); 

K v = k/(v F A), 

1 rl 

Ki(k v ,\) = — dx\J\ — 
IX Jo 



X 2 



\(v 2 + x 2 ) 



V 2 — X 2 ,0/2 2\ 4AXZ/ 7T K v 

2\u In A (k v + x ) - -(- - arctan — ; 

is z + x A v A + x A 1 \x\ 



K 2 (k v ,X) = — / 

7T JO 



dx 



VT^ 2 

k! + x 2 ' 



v = K v - -. (4.15) 

Below we mainly use k v instead of k. We find that no terms proportional to sgn(a;) ap- 
pear in the above expansion, and so the renormalization constant of k, R K (t), equals to 
unity. Therefore the renormalization of k comes only from the wave function renormal- 
ization R(t). From ( 4.14j ), we obtain the renormalization constants for infinitesimal RG 



transformation as follows 

Rn{dt) = 1, 

R(dt) = 

= 1 + adtKi(n v , A), 

R VF {dt) = 1 + S W + S P 

= l-adtK 2 {K v ,X). (4.16) 

In Eq. (|4.15|) the dominant contributions to the integrals for K\ and K 2 come from the 



region x ~ 0, so we can approximate as yl — x 2 ~ 1 in the integrand. Then the x 
integrations can be performed readily as 

2 11 1 
Ki(k v ,X) ~ —arctan — — (in A 2 + ln(K 2 + 1) — 2i/arctan — 

7T K v 7rA(l + V 2 ) V ' K v 

2 1 

K 2 (k v ) ~ —arctan — . (4-17) 

7T K n 
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4.2 Renormalization of the gauge-field part 



By integrating the high-momentum modes, we obtain the one-loop correction to the action 
Sa, so the renormalized action Sa becomes 

p A(l-dt) r o n | e | 



A^veiq 2 -' + A 2 - 6 A^) + U(q, e) A(q, e), (4.18) 
VFq J 



Sa = d<L dqA(—q, — e) 

where the vacuum polarization II (q, e) reads explicitly as 

n(g,e) =9 2 4E(^)7^ / A *G (a;|) + ^ + e)Go(a;^). (4.19) 
~r v o 7 j JA<i-dt) — 



g / j ja(i-m) 

After doing the integration in the radial direction of g, IT(g, e) is expressed as 
g 2 v F Adt M-^fe a x q\2 rv F A 



n(g,e) 



2tt 2 



V v <? 7 ^-«fA |sm0(y)| 



(4.20) 



where 



3 3 

P(e,E q ,y, K ) = e(e)X)P i (e, J E?„ J /,«) + e(-e)X;P i (-e,- J E?„-y,«), (4.21) 

i=l i=l 



with 



P!(e,E q ,y, re) 
P 2 (e,E q ,y, re) 
P3(e,E q ,y, re) 



/ 

J—oo 

/ rfo; 

J-E 



i{uj + e) — ik 


-y- 


- E q iuj — in - 




1 




1 




+ e) + m - 


-y- 


IUJ — IK ~ 


K 


1 




1 




"i(c<j + e) -Hre 


-y- 


- EqlUJ + IK - 


-E q 



(4.22) 



We have changed the integration variable from to y as before. The above integrations 
can be evaluated as follows; 



Pi(e,S„j/,K) = 



P2(e,E q ,y,K) = 



2ir(ie - E q 



1 re 2 + £ 2 

2 n (e + re) 2 + £ 2 



«sgn(E)(| -arctan^; 



-Msgn(y)( arctan- 

2 1 2/ 1 



2vr(i(e + 2re) - 



1 (e + .) 2 + £ 2 _ ln re 2 + £ 2 
2 l re 2 + £ 2 (e + /t) 2 + i/ 2; 



—isgn(E) (arctan — ^j— arctan — — ■) 

\E\ \E\ 

+isgn(y) (arctan — — arctan ) 

\y\ \y\ 
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2n(ie - E q ) 



— In - 

2 k + y 



2 



— isgn(E)(— — arctan— j—) + zsgn(?/)( — — arctan-—) , (4.23) 
2 \E\ 2 \y\ - 1 



where E = E q + y. We are interested in the low-energy modes, so we expand these results 
in powers of q and E q . Detailed calculations are given in Appendix B. Explicitly, for 
\ 6 \/( v fq) << 1 @, (J we obtain 

n(5 , £) = ^ E J"L_f t 4 «• 



where 



2tt 2 V 27r(ze - E q ) Jo ^T^? «2 + 3:2 

. 7rA;p lei 2 1 

VBAaat— arctan — , (4.24) 

A v F q 7r «;„ 



K a = (^«) 2 (4.25) 



and we set \/l — x 2 ~ 1 in the x-integration above as before. This result has a form of 
the dissipative (A) term. We neglected a constant term, which is a mass renormalization 
of gauge field, because this term is canceled by the contribution from the four-point 
interaction as explained in the introduction. 

By substituting these results into Eq. (|4.18 ), we get the infinitesimal renormalization 



constant for the coefficient of the dissipative term, 

R x {dt)X = A + adtK 3 {K v ), (4.26) 

where 

K 3 (k v ) = — ^ arctan — . (4.27) 
A K v 

4.3 Vertex correction 

The renormalization of the vertex part is similarly obtained by the momentum-shell in- 
tegration, and the action Si nt .3 is expressed as 

rA(l-A) j j ^ ea x q 

q 

x A(g, e)tp(a;p + q,u + e)ip(a;p,u), (4.28) 



Sint.3 = igv F J dude J dpdq £ a 6a * g (l + T 3 (a; p, u, p + q, u + e 
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where the vertex correction T 3 is given by 



r ( ' >\ a '9 v F(e a x q\2 

r 3 (a;p,o;,p,w) = / dq ——[ 

JM-di) — 1W v a' 



/A(l-cft) — l{q') v <?' 
x U(E p + E q ,,uj;E p , + E q ,,uj';F(q')) } (4.29) 



and 



U(E,w,E',w'\ F) 
= /'^k'G'o(a; J E,^ + e')G' (a; J E / ,cu / + e , )^ 

| C | "( 

= / <k'(G (a;£,cc> + e ') - G ( a; #V + e'))r^ 

J £ ~ | 



+ F 

~+F 



1 



(4.30) 



£(c</ — oj) + m(sgn(a/ + e') — sgn(u; + e')) — (E f — E) 
The corrected three-point coupling has a momentum dependence. However, we take the 
"on-shell" limit since we are not interested in irrelevant higher-order terms of the Taylor- 
expansion in powers of momenta and energies. The on-shell limit of T 3 should be taken 
as § 

lim hp r 3 (a;p,w,pV) = 1™ F 1 F (S(0, - S(0, ^) + 0(E 2 p , E 2 p ,)) 

= S P , (4.31) 

where S p is given by Eq. ( |4. 14| ) . Therefore we obtain the renormalization constant R g (dt), 



R g {dt) = 1 + S P 

= R(dt)R VF (dt). (4.32) 

This relationship between R g = RR Vp is nothing but the Ward-Takahashi identity which 
is discussed in Ref.[l5i. 



5 WRG equations and the RG flows 

In the previous section, we obtained infinitesimal forms of renormalization constants, 
Eqs.( |4.16| ), ( |4.26| ) and ( f4.32| ). These results and the scaling laws fl3.4|) give the differential 



14 



equations for the running parameters, g(t),vp(t), X(t) and n v (t) as follows, 



dg{t) 
dt 


2 - 


2 SW ' 


dvp{t) 
dt 


(« 


, s(*)V(i) v , 4lV „ ,^ 


d\(t) 
dt 


( 2 


_ 6 _ 9 WVW A - (t))A(i) + 9 WV(0 /f3(t)i 


dn v {t) 
dt 


(i 




i = 1, 2 and 


3) 


means X(t)) given by (|4.17) and (|07| 



(5.1) 



In these equations, in place of the gauge coupling g itself, one can regard the following 
combination 7 as an effective expansion parameter in the perturbative calculation of the 
present theory, 

7 s £S£!lW (52) 

In Ref.|| we have introduced a similar combination ( we called it a = g{t) 2 v F{t)(2n 2 v 
there), where k = and so K 2 (t) = 1. This 7 is a generalization of it in the case of k > 0. 
The reason why the extra factor K 2 (t) appears is physically obvious from (|1 . 5|) and Fig.l. 
That is, the density of states around the Fermi surface decreases by the existence of the 
K-term as N(Ek) ~ 1/2 ± (arctan k~ 1 )/-k ~ (1 ± K 2 )/2, and so the interaction between 
the gauge field and fermions is weakened by the factor K2(t). 
Then we obtain the following coupled differential equations, 
dj(t) 



dt 



1 _ b _ Ll{ t) - (1 - L 1 (t)L 2 (t)) 7 (t) 7 (t), (5.3) 



= {i-l-lit))v F {tl (5.4) 

^ = (l-L 2 (tHt)) Kv (t), (5.5) 

= (2-6- 7 (t))A(t) + ^7(t), (5.6) 
where -^i(t) and i^2(^) mean the following L%(K v (t) , X(t)) and L 2 (K,(t), A(t)), respectively; 



1 ^ (k^ + 1) arctan k v 1 ' 
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For latter discussions on the IR behavior of the parameters, it is somewhat useful to refer 
to the system with k v — and K 2 — 1. Therefore we rewrite below these differential 
equations in terms of the combination, 

-M - _ifiL = ^M! ( (5 .8) 



l-&-AT 2 (t)a(t))a(t), (5.9) 

1 - {K x (t) - K 2 {t))a{t)) Kv {t), (5.10) 

2 - 6 - K 2 (t)a(t))A(t) + K 3 (t)a(t). (5.11) 



as follows; 

da;(i) 

dn v {t) 
dt 

d\(t) 
dt 

The equations (|5.3|), ( |5.5| ) and ( |5.6| ) are a set of closed differential equations with 
respect to the coupling constants, j(t),K v (t) and A(i), and they do not depend on the 
scaling parameter £. The behavior of Vp(t) is calculated from these quantities. Then, in 
the following discussions, we shall consider RG flows of only these couplings. 

As stated, we want to know whether the present theory including the K-term has a 
nontrivial fixed point or not. In the case k v = 0, this theory has the nontrivial IR fixed 
point in the parameter region b < 1 such as 7(00) = 7q = 1 — b > 0. The existence of 
this fixed point can be seen in Eq. ( |5.9|) by simply setting i^ 2 (^ = 0) = 1, where k v = 
is certainly a solution. We shall discuss how the existence of the K-term influences the 
stability of this fixed point. 

A formal solution to the differential equations ( |5.3|) , ( |5.5|) and ( |5.6|) is obtained as 

= e^K 2 (t) 
7U j t Q dse^- b >K 2 {s) + a^' 

K v (t) = K v0 exp(t- dsL 2 (s)i(s)), 

Jo 

X(t) = A e (2 - b) ^o^) + ^. (5.12) 

One can determine the qualitative IR behavior of the parameters by the following argu- 
ments; 

(i) First, let us assume k v has a finite constant value in the IR limit, i.e., lim^oo K v (t) = 
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k* > 0. Then, from Eq.( |5.9| ), a(t) has the following fixed points a* 



l-b 



(6<1) 



a* = { K2(K) v . (5.13) 

(1 < 6 < 2) 



Then X(t) becomes large according to ( 5.11|) as 



A ~ (2 - b - K 2 a*)X + K 3 a* 

f A + K 3 a* > (b < 1) 

= (5-14) 
I (2 - 6) A + K 3 a* > (1 < 6 < 2), 

where A = d\(t)/dt. Since A becomes large, the factor K\ — behaves like, 

In A 2 

K x -K 2 ~ ^^0. (5.15) 

7rA 

This fact and Eq. (|5.10|) indicate that k v increases exponentially. This is an obvious 
contradiction to the starting assumption. This implies that k v cannot have a finite value 
in the IR limit. 

(ii) From (i) we conclude that the parameter k v tends to large in the IR limit. Thus 
L\ goes to unity and L 2 is given as follows, 

1 



K~ 2 X 2 + (x — l) 2 



[xlnx + 1 — x), (5.16) 



where x = k v X. We note that the value of L 2 is restricted to — 1 < L 2 < 0. The coupled 
equations are approximated as follows, 

7 ^ (-6+ (L 2 - 1)7)7, 
k v ~ (1 - L 2 j)k v , 

A ^ (2-6- 7 )A + ^7, (5.17) 

where we have omitted the argument t in j(t) etc. In (ii-1,2,3) below we argue that 
j(t) ^ in the IR limit. 

(ii-1) Let us assume that j(t) has a nontrivial fixed point 7* = j^zi > 0- Since ^(t) 
should be positive, L 2 should become larger than unity. However we found that L 2 < 
in (ii) above, and so this cannot happen. 
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(ii-2) Let us assume that ~f(t) tends to infinity. Then, A(i) satisfies the equation 
A ~ — 7A + (7r/ci?/A)7, and we obtain \(t) as 

roc 7lkr? 

X(t)=exp(-J o dsj(s)) + ^-. (5.18) 

From the behavior of K v (t) and X(t) at large t, k v — > 00 and A — > nkp/A (constant), we 
find that L 2 ~ —Inx/x — > 0. This gives 7 ~ (—6 — 7)7 and 7(t) (7(£) should be positive) 
tends to vanishingly small. This contradicts the starting assumption. Therefore, j(t) 
cannot tend to large in the IR limit. 

(ii-3) After all, the only possible situation is that 7* = 0, i.e., IR limit is the trivial 
fixed point. In this case, k v and A become large (for < b < 2) according to the original 
canonical scaling such as 

A ~ (2-6)A. (5.19) 

Since we have obtained the qualitative behavior of the parameters in the IR limit, let 
us turn to the numerical calculations of the RG flows. The result shows not only that the 
above qualitative argument is correct but also that some interesting crossover phenomenon 
appears. We discuss the two cases (1) < b < 1 and (2) 1 < b < 2 separately. 

(1) The case < b < 1. 

When k v = 0, as mentioned, 7(i) has the nontrivial IR fixed point 7q = 1 — b > 0, 
and at this fixed point the system behaves as a non- Fermi liquid^, |J. On the line of 
k v = 0, the RG flows converges to this point. The RG flow whose initial value of n is 
very small, k v ~ 0, first approaches the point Po{k v — 0,7 = 7q > 0). Then k v increases 
and the flow goes away from the point Pq, finally approaching the trivial fixed point 
Pi{ftv = 00,7 = 0). Therefore, the point Pq is an unstable fixed point of this theory and 
the IR limit of the system is controlled by the point Pi which describes a Fermi liquid with 
impurity effects. Such behavior is shown in Fig. 4. and Fig.5. In Fig. 4 we plot the flows 
close to 7(i) = 1 — 6. In Fig.5 we show the global behavior of RG flows in a wider region 
of k v . In Fig. 6 we show the wave-function renormalization constant R(t), together with 
the coupling constant j(t). As long as j(t) is staying in the vicinity of the unstable fixed 
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point P , R(t) increases exponentially ~ exp(2r/t), where rj is the anormalous dimension 
of fermion field (for k v = 0). As discussed in Ref. 0, this behavior of R(t) implies that 
the system in these energy scales is controlled by the unstable fixed point Pq and it has 
non-Fermi-liquid-like behavior. This result will be confirmed by the discussion on the 
fermion propagator which is given in Sec. 6. As RG flow is going away from this unstable 
fixed point, approaching the trivial fixed point Pi, 7 decreases quickly and R(t) becomes 
constant, as shown in Fig. 6. 
(2) 1 < b < 2. 

The RG flows of this situation are shown in Fig. 7. They reach the point Pi in the IR 
limit. In both and A increase exponentially. 

6 Low-energy behavior of fermion propagator 

In this section, let us consider the behavior of full fermion propagator under the RG 
transformation with the generalized scaling assignment (parameterized by £): 

u = e^uj, p = e'p, 
ip{a;p,u>) = e^ +1)t {R(t)}^(a;p,uj). (6.1) 

We begin with the definition, 

(if)(a;p,u)if)(a';p',u/)) = -5 a>a >5(p - p')6(u - u') 

x G(p, uj; vp, 7, A, k; A). (6.2) 

By substituting the last line of Eq. ( |6.1|) into the LHS of Eq. (|6.2| ) and then expressing it 
by the propagator at p, we obtain 

Ma;p,^(a';pV)> = e 2 ^ {R(t)y l ^{a; p, S)^(a'; p', a>')> 

= -e^+^iR^rXJip - ff)S(Q - Co') 

xG(p, u- v F {t)n{t), A(t), A). (6.3) 

From (6.2) and (|6.3|) we have 

G{p ) u- Pl -K) = ^ t {R{t)Y l G{p,Co- Pl {t)-K) (6.4) 
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where Pi(t)'s denote the running parameters VF(t),y(t), X(t) and n(t). This equation can 
be used to evaluate the full fermion propagator simply by substituting the solutions of 
the RG equations to the RHS of Eq.(|6.4|). 

Explicitly, the RHS of Q6.4| ) is written as follows; 

e^ t R(t)G- 1 (p,Co;p l (ty,A) 
~ e~^R(i) iuo + m(t)sgn(tl>) — v F (t)e a ■ p — £/(0, u; Pi(t); A) 
~ iR(t)uj + iR(t)R K (t) Ksgn(u) — R(t)R VF (t)v F e a ■ p 

—e~®R(t) £/ (0, u>; pi (t) ; A) , (6.5) 

where £/(]5 = 0, Co; Pi(t); A) is the fermion self-energy which is obtained by integrating out 
all the modes of fields (i.e., with momenta smaller than the cut off A). We have neglected 
its momentum- dependent part which is less dominant than v F e a ■ P at low energies. We 
evaluate the self-energy Sy(0,u;; pf, A) at one-loop level as follows; 



£/(0, A) 



wbA 



e a x q 



i(u + e) + wsgnfa; + e) - v F (t)e a ■ q ( £\ 2 , A(t)|g| 

; 9 2 (t)v F (t) f ,„ f A sgnp + e) 
' (27r) 2 w B A 



" (*) 



+ 



v F (t)g 



\u + e + nsgn(u + e) 
v F (t)q 



(a; + e) + /tsgn(u; + e) 
v F (t)q 



(6.6) 



where we have performed the angle integration, and used the fact sgn(u>+e-|-fi;sgn(u;-|-e)) = 
sgn(u) + e). In the e-integration above, the region of e ~ gives dominant contribution, 
so we treat e as small and restrict the integration region to < e < v F (t)q. We also treat 
Co small because we are interested in the low-energy behavior of the propagator. Below, 
we discuss the three cases, (I) < 6 < 1, (II) 6 = 1, (III) 1 < 6 < 2 separately. 
(I) < 6 < 1 

In the previous section, we showed the RG flows of the parameters j(t), K v (t) and X(t). 
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In the small K v (t) region, the expansion parameter ^(t) approaches the saddle point P 
with 7q = 1 — b. As K v (t) tends to large, the parameter ^(t) leaves away from P and 
approaches the trivial fixed point P±. We study the low-energy behavior of the propagator 
in two regions; (i) region near the saddle point Pq where K v (t) is small, and (ii) region 
near the trivial fixed points P\ where K v (t) is large, 
(i) Small K v (t) 

For small K v {t), the self energy is evaluated as 



E,(0, A) 



2AK 2 (t) J 



de / dq 
Jo 




(i\ 2 - b + mm 

\AJ ^ v F (t)q 



1 



sgn(w + e) 1 



u + e| + n(t) 



v F (t)q 



) 



+ 0(iU 2 ,K v 2 (t)) 



(6.7) 



where 




l K 2 (t)v F (t)AJo € Jo q q 3 - b + X v (t)e 




(6.8) 



and 



X v (t) 



X(t) 



v F {t)A 




(6.9) 
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The upper limit of integral region with respect of e in S 3 is vp(t)A then the argument 
contained in S3 is X(t) and is not X v (t). To obtain an explicit form of the propagator, we 
need to evaluate these if's. 

Now we consider the case where t increases so that X(t) becomes large, while j(t) is 
still staying near the saddle point and the argument c = A„(i)|u)| of Si, S 2 is still small. 
Then, we evaluate the above integrals for small c 

H u (c) ~ 

H 2b (c) ~ ^l_(l- c §% (6.10) 

On the other hand, the argument X(t) of if 2ft i n ^3 is large, and then for large c, 

H 2b (c) ~c^(lnc + 3-6). (6.11) 

We have used the assumption of small A v (t)|a)| above just to evaluate the integrals Hi b 
and H 2b . The discussion later on shows that the Green function obtained below by using 
( |6.10| ) is not only a solution to the RG equation at low energies for small A„(t)|u)| but also 
remains to be an approximate solution even for large A„(t)|u)|. 

By using ( pMOj ) and (|6.11|) , we obtain the Green function as follows; 



e-^R(t)G-\p,u; Pl (ty,A) 
= iR(t)u + ie-^R(t)K(t)s gn(cj) — R(t)R Vp (t)vFe a ■ p 



+ K^-r 1 i; + 4 ( " lA(i)+3 - 6) ))^ ( " 2) 

In the situation where X(t) becomes large, while 7(i) is still staying near the saddle point 
and k v is small value, we obtain the folloing equation from (|4.17| ) such as 



K L ~ K 2 = 1-^ + 0(k 2 v ). (6.13) 

7T 

From this relation and R(t) , R VF (t) in Eq. (f4.16|) , we obtain 

R(t) ~ e/o^W 
R(t)R VF (t) ~ 1. (6.14) 
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The formal solution of j(t) which is given by fl5.12| ) and the above equation give the 
following relation, 

R(t)^{t) = « e (1 - 6)t K 2 (t), (6.15) 

where «o is the initial value of a(t) of ( |5.8|) . From (|5.1|) , we also find 

v F (t) ~ e (5_1) Vo 
n v (t) ~ e l K vG < 1 

\ v (t) ~ A, e (3 - 6 -^. (6.16) 
By using these equations, i?(t) is obtained as follows; 

R(t) = 1 + a J* dse {1 - b)s K 2 {s) 

= 1 + Ti( e(1 - s, '- 1 )-^f)(^"- 1 )- <« 7 > 

We substitute these relations, which depend on £, to Eq.( |6.12j ), and obtain the Green 
function as 



G-\p, u- Pl - A) = e-^R(t)G-\p, Co; Pl (t); A) 
~ ia sgn(u;) 



/ 1-6 2-6 
/ \ 3-6 „ \ 3-6 



\v I i 2 Kj V \d 

a; 3-6 Icul 3 - 6 



1 - b 1 1 2-6 



+ ^ ( i - + ^(2 - 6) / + lKS & n ( u ) ~ Vp6a ' p + A(j 

Gleading) _1 +OH, (6.18) 



where 



AG- 1 = iua [^^bt^ - 1) - ^(ln A + t + 3 - 6)V (6.19) 

and we have omitted the suffix for initial values of a, A, X v , k v and vp. 

Here we comment on the last term AG -1 which depends on the RG scale t. The LHS of 
Eq. ( |6.18D is independent of the scale parameter t, so the representation on the RHS should 



not have terms depending on t if we calculated the fermion self-energy exactly. However 
one cannot do the exact calculation, and we evaluated the self-energy perturbatively at 
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one loop level assuming the RG parameters lie in the region where 7(1) < 1, If this 
region corresponds to the infrared limit (like the case k v ^> 1 and 1 < b < 2 which we will 
discuss later), the t-dependent terms should disappear in the infrared limit (t — > 00) with 
exponentially dumping factors, for example see Eq.( |6.25[ ) and Eq. fl6.30p . In that case, the 



one- loop calculation is almost exact. Right now, we are considering the situation with 
some finite t, t ~ £*, in which j(t) is staying near the saddle point 7q = 1 — b. We 
assume 1 — b <C 1 so that the perturbative calculation assuming small 7 is meaningful. 
However, the above AG -1 does not disappear even in the limit t — > t* with the t dependent 
factor such as e^ 2 ~ b ^. Of course, we find that the leading terms of Eq. (|6.18|) , Grading; are 
independent of t and are in fact a solution to Eq. ( |6.4|) up to the first order of K v (t) 
and at low energies. This fact suggests that the AG -1 term shall be canceled out with 
higher-loop corrections to the fermion self-energy. Actually, the two-loop corrections are 
known to contain terms of the order (7q) 2 /(1 — b) = 7q, which are of the same order in 
7q = 1 — b (<^C 1) as the one-loop correction. 

From ( |6.18| ), we find that the dominant terms of the Green function (Gi ea dmg) are 

2 1 

0{ujz-b) and 0{k v uj^) in the present parameter region < b < 1. Therefore, at interme- 
diate energy scales, the Green function has the behavior of a non-Fermi liquid; a branch 
cut rather than a pole in uj appears as in the Luttinger liquid in one dimension||]. 
(ii) Large K v (t) 



In this case, the self-energy of (6J;) is evaluated as 



■ l{t) , ^ /V f A ^ 1 v F (t)q 

— —i—. r^sgn [uj + e de dq ^—7 — 

2AK 2 (t) B K J J Jo (l\ 2 ~ b , MM. 2/c(t) 

VAj ^ v F (t)q 



3 



^MiU ^ 1 '^'^'^'' (6 ' 20) 



where 



H 3b (c) = — I dyy^ ln(l + y). (6.21) 



For small c, is evaluated as 



H 3b (c)~~(l-c^). (6.22) 
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Then we obtain G 1 (p,u; pf, A) as 

e-« t J R(t)G- 1 (p,w; i o i (t);A) 

~ iR(t)uo + msgn(w) - w F e a • p + ie~ ft ^ \; 7 ® a) f (A„(t)p|)^ - l) (6.23) 

where we have used K 2 (t) ~ 2/(7TK v (t)). Now k is large enough, so behaves as 

2ao /■* 1 

(is — — 



JKt) = i + 

7T JO 



K V {S) 

~ 1 + 70(1-6-'), (6.24) 
where 70 = 2ao/(vr/t 1 , )- Thus, in this case, the Green function is obtained as 
e-^R{t)G~\p^- Pl {t)-K) 

~ i/csgn(w) - fire a ■ p + i (1 + 7) a; — -^-Xv~ b sgn(u;)\u;\3=b + 0(e~*) (6.25) 

where we have omitted the suffix as before. 

In the low-energy limit, we find that the dominant term is linear in u, so the Green 
function has the Fermi-liquid behavior. This result of course agrees with the conclusion 
in Sec. 5 obtained based on the RG flow. 

(II) b = 1 

In this case, 7(t) goes to the trivial fixed point P\. 
(i) Small K v (t) 

For small K v (t), we obtain the self energy by setting b = 1 in Eq. fl6.8Q . The functions 
H\i and if 21 are easily obtained as 

H n {c) = I((I + l)ln(l + c )-]nc) 

H 2 i(c) = -^=In(l + c) + 2arctan^. (6.26) 
V c V c 



From these equations and the t-dependences of parameters in ( |6.15|) , ( |6.16|) and ( |6.17|) 
with b = 1, we obtain the Green function as follows; 

~ iKSgn(u) — VFe a ■ p + i(l H — + — lo In (A„|a;|) + 0(t) (6.27) 

7T 2 v ' 
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where the suffix for initial values are omitted as before. The last term 0(t) is the t- 
dependent part which should be canceled out with higher-order corrections as explained 
above. 

This expression contains ujln{u) term, leading to the weight of quasiparticles that 
vanishes logarithmically at low energies. This behavior is like that of the marginal-Fermi 
liquid @. 

(ii) Large n v (t) 

As t becomes large, n v (t) becomes large too. Then we evaluate the self energy as 

£(0,u>; A (t);A) ~ -isgn(a)) f .. Mt)\Q\)iH al {K(t)\a\). (6.28) 

Here H31 is evaluated as follows 

c*H 31 (c) = - (ln(l + c) + 2c(l - v^arctan ^=) j . (6.29) 



By using ( p.lC ) and ( 6.24 ) with b = 1, we obtain the Green function as 



e-^R{t)G~\p, u>; Pi {t)- A) ~ msgn(cu) - v F e a ■ p + i(l + 7)0; + 0(e - *). (6.30) 

The logarithmic term has disappeared but only the term linear in u survives (together 
with the signature term). This behavior of the Green function is that of the ordinary 
Fermi liquid. 

(Ill) 1< b < 2 

In this parameter region, the Green function has the same form with (|6.18 ) for K v (t) 



1 and (|6.25 ) for n v (t) ^> 1. But the parameter b > 1 here, and so the dominant term 



of the Green function is linear in uj at low energies. Therefore, we find that the Green 
function has the Fermi-liquid behavior in both cases of K v (t) <C 1 and K v (t) ^> 1. 



7 Concluding Remarks 

In this paper, we have given detailed studies on low-energy behavior of the system of 
nonrelativistic fermions interacting with a dissipative gauge field by using the WRG. 
Especially, we investigated the stability problem of the nontrivial IR fixed point in the 
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clean system against impurities. The K-term in the fermion propagator, which is generated 
by the interaction between fermions and impurities, is introduced for this purpose. We 
showed that this term makes the non-Fermi-liquid fixed point unstable, and in the IR 
limit the effective gauge-coupling constant tends to vanish. However in intermediate 
energy scales, the system is controlled by the unstable fixed point as a saddle point, 
exhibitting non-Fermi-liquid behavior for some parameter regions of b and k. This result 
is confirmed by obtaining the fermion propagator by solving the RG equation and also by 
calculating the wave function renormalization constant. 

It is also quite interesting to investigate finite-temperature effects in the present gauge- 
fermion system by the RG analysis; in particlular, the same stability problem. This 
subject is under study and we plan to publish the result as a paper subsequent to the 



present onellO]. 
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A Taylor series of the self energy 



In this Appendix, we make the Taylor expansion of the fermion self energy () to obtain 
the result (|4.14j) . It is assumed that Zq = u + E p takes a small value and the self energy is 
expanded around y = E q . After the y integration, only even functions of y in remain, 
and then we neglect odd functions of y in the expansion of 

t / n \ 1 / 1 1 (u; + l) 2 . . . ,7T K . 

h(u,E,K) = -^ TI (2 ln ^^ + 2Sgn(E)( 2- arCtan M 
= c (K,y) + d(K,y)z + c 2 (K,y)iE p , 

I 2 (u,E,k) = — — -(-(ln^ _ F9 + ln|io + l| 2 ) 



Co" + K K 

zsgn(i?)(arctan — — — arctan —— 

\ E \ \ E \ 
di(K,y)z + d 2 (n,y)iE, 



p- 



where 



I 3 (co,E,k) = -^(^(ln((w + /€) 2 + £ 2 )+*sgn(£)(^- arctan 

= 9q(k,v) +9i(K,y)z + g 2 (K,y)iE p . (A.l) 



Co(K,y) = ^((1 - k) ln(K 2 + y 2 ) - 2|y|(^ - arctan-^; 
ci(k,j/) = — («- 1 + Z>(k, y) 



2 



and 



a. 

ca(«,y) = - — ((l-«)(-^-5-l) + ^-^)+S(K,y), 
a_ v «r + k z + y z/ 

di(K,y) = -d 2 (n,y) = K 

k z + y A 

0o(«,y) = -Co(«,y), 

ci_ * K -\- y ' 
g 2 (K,y) = E{K,y), (A.2) 



D(K,y) = - — j (-c- (k, - ln(/« 2 + y 2 ) + 4(1 - K.)\y\Q - arctan ■ , ,. 



2(1 — K) /7T K \ 

E ^ = -w^YU- arctan u)' 



a-(«,y) v 2 ' |y|- 
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a-(K,y) = (1 - «) +y , 

c.(K,y) = (l-K) 2 -y 2 . (A.3) 
The self energy is expanden as follows; 

S(a; p, uj) = S + S^zcj + S p £ p . (A.4) 

where 

S = -iadtF A f dy\ sin (p(y)\sgn(cu) (c (k, y) + g (~K,y)) 

J — A 

= o, 

A 



adtF A / dy\ean<j>{y)\(ci(K,y) + d 1 (K,y)+g 1 (K,y) y j 

J — A 

1 f 2{K - 1) - C -K y \ H- 2 + V 2 ) + Alyl{ }~ K \- - arctan * ] 



a_(K,y)V a-(«,2/) o-(« s y) 2 M' 

S p = / p — S w . (A.5) 



and 



Jp = adtF A f dy\ sin 4>(y)\(c 2 (K,y) + d 2 (n,y) +g 2 (n,y) 

J — A 



-A 

* -2k 



= adtF A / efy|sin0(y)| _ 2 , (A. 6) 

./ —A k \ y 

y and 7c are defined in Eq. ( |4.7|) such as y = VpA cos /Fa and 7c = k/F a . 



B Vacuum porlarization 



In this Appendix, the vacuum porlarization ( |4.2G|) is calculated to give ( 4.24j ) . In the y 



integration, only even functions of y in Pi survive, and odd functions of y are omitted. 
Then Pj's are obtained as follows, 

P2(e,E q ,y,K) = (- ln ^- - - -j—) iE qj^ 

P 3 (e,E q ,y,K) = (l\n^- 6 -±^iE q ), (B.l) 



1mA 








2 L e 




-InA 


k e + k 












2 L e 
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where 



L = k 2 + y 2 , L e = (e + nf + y 2 . (B.2) 

We assume that the coefficient of the impurity term k is larger than the momentum and 
energy, q and e, as ij/k C 1 and e/k C 1. In this situation, we find that Pi and P 3 are 
order 0(1/ k) and P 2 is order 0(1/k 2 ), and then we can neglect the P 2 term. We also set 
ln(L/L e ) ~ because L ~ L e . After all we obtain the approximate form of Pi as 

£l\(e,E 9 ,y,K) - FT (y + C + /6x 



i= i 27r(ie-P g ) L L e 

Eg « 



7r(ie — P g ) L 

The vacuum porlarization part is given as follows 



2tt 2 V 27r(ie - £,) Jo ^/l^x 2 k 2 v + x 2 
. , kp r n dd) vfQ cos d> sin 2 6 f 

~ VB^Oldt— I ; / 

2A J —it 2ti ie — VfQ cos m Jo 



2ir ie — vpq cos dj Jo y/1 — x 2 n 2 + x 2 
Icf \e\ f 1 , 4 



~ const + vshadt-^-^- [ dx = v 
2A v F q Jo V 1 - x 2 ni +■ 



where 



and use the relation X^e a (f^) — S-w^- ^ n the x integration, we set y/1 — x 2 ~ 1. 



(B.3) 



7r/c F e 2 1 
const + Vs-h-ddt— arctan — , (B.4) 



/C = (^) 2 , (B.5) 
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Figure Captions 

Fig.l 

Momentum distribution of the fermion density. The existence of the K-term smears the 
step-function-like behavior at k — 0. The energy scale corresponding to the momentum 
cut off A is also shown. 

Fig.2 

— * 

Partition of the two-dimensional space of fermion momentum K . The big arc C—C shows 
a part of Fermi circle of radius K F around the origin O. Theare are TV = 2irK F /(2A) 
Fermi wave vectors, K Fa (a = 1, ...,7V), being separated each other by equal distances 
2A along the circle. Here A is the cutoff momentum. Under the renormalization-group 
transformations, we integrate gradually the modes with the momenta that sit "near" the 
Fermi circle within the distance ~ A. To be explicit, we integrate the momenta wich lo- 
cate within one of the TV small circles of radius A, each circle being around the end point 

— * 

of K Fa as illustrated. One may enlarge the region of relevant momenta, e.g., to those 
that locate within the union of small circles of radius 2A, to avoid the momenta which sit 
very near to the Fermi circle but never be integrated. However, the main conclusions of 
the text do not depend on such detail. 

Fig.3 

The two-dimensional space of momentum k = (k x , k y ). k = K — K Fa is a deviation of 
fermion momentum K near the Fermi circle measured from the nearest Fermi momentum 
K Fa . In the step of renormalization-group transformation from t to t + dt, the modes 
with such momenta that sit in the shaded area are to be integrated. This area is nothing 
but a shell surrounded by two circles of radii Ae~* and Ae~( t+dt \ 

Fig.4 

Structure of RG flows of K v (t) and j(t) near the unstable fixed point Pq (7*, k v ) = (1—6, 0). 
Values of the parameters are b = 0.5, k F — 1, A = 1 and A = 0.1. For sufficiently small 
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initial value of k v0 , the RG flows approach to the fixed point P . 



Fig.5 

Global structure of the RG flows. Same as Fig.4, but the scale of k v is enlarged. The 
effective coupling constant j(t) tends to vanish in the IR limit. 

Fig.6 

The solid line is the effective coupling constant 7(t), and the dashed line is the renormal- 
ization constant In R(t). As long as the effective coupling constant stays near the unstable 
fixed point Pq, the renormalization constant grows up exponentially. 



Fig.7 

RG flows for 1 < b < 2. K v (t) increase exponentially and 7(2) tends to vanish in the IR 
limit. 
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